S-wave superconductivity with orbital dependent sign change in the checkerboard 

models of iron-based superconductors 
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We study three difTerent multi-orbital models for iron-based superconductors (iron-SCs) in the 
solvable limit of weakly coupled square plaquettes. The strongest superconducting (SC) pairing is 
in the Aig s-wave channel and its development is correlated with the emergence of the next-nearest- 
neighbour antiferromagnetism (NNN-AFM). For the models with more than three orbitals, this 
study suggests that the signs of the intra-orbital pairing order parameters of the dxy and the dxz (or 
dyz) orbitals must be opposite. Such sign difference stems from the intrinsic symmetry properties of 
inter-orbital hoppings and might, ultimately, lead to the sign-change of the SC orders between the 
hole Fermi pockets at the F point and produce anisotropic or even gapless SC gaps in the electron 
Fermi pockets around the M point in reciprocal space, as restoring back to the homogeneous limit. 



Introduction- Since the iron-based superconductors 
were discovered two years ago 5]| , the relation between the 
new superconductors and the high- Tc cuprates has been 
a central focus of researches. It is highly debated that 
whether the new superconductors belong to the same cat- 
egory of strongly correlated electron systems in which the 
cuprates are believed to be. Models based on both strong 
coupling and weak coupling approaches have 

been applied to understand the properties of the new ma- 
terials and their relation to the curpates. 

Strong correlation in an electron system can be re- 
flected in its strong "locality" of physical properties 
caused by short-range interactions. For instance, both 
magnetism and superconductivity in the cuprates could 
be attributed to such "local" physics. The magnetism in 
the parent cuprate compounds is well described by the 
Heisenberg model with the nearest-neighbor (NN) anti- 
ferromagnetic (AFM) exchange couplings. The d-wave 
superconductivity with symmetry form factor, cos kx — 
cosfc,,, in reciprocal space [17| corresponds to short-range 



superconducting (SC) pairings between just the NN cop- 
per sites. Theoretically, the locality (short coherence 
length) already allows us to understand many essential 
physics from small clusters and to give new insights for 
obtaining low-energy effective models 181-121 1 . 

Remarkably, recent experimental results in iron-based 
superconductors suggest that both magnetism and su- 
perconductivity display an excellent locality behavior as 
well. Neutron scattering experiments demonstrate that 
the magnetism in the parent compounds of iron-SCs can 
be again described well by the Heisenberg model with the 
NN and the NNN-antiferromagnetic exchange couplings 
between iron spins (22| - [23 |. In addition, angle resolved 
photo emission spectroscopy (ARPES) experiments show 
that the SC gaps can be fitted to a single cos kx cos ky 
functional form in reciprocal space j25l427l |. Regarding 



the fact that iron-SCs are complicated multi-orbital sys- 
tems, these experimental results are compelling evidence 
supporting that a local interaction controls the pairing 
channels since the function form in the real space simply 
corresponds to SC pairings between the NNN iron sites. 
They motivate us to ask whether the physics of iron-SCs 
can also be understood from checkerboard-like models 
as studied in Ref. pjj, where the local physics (within 
a plaquette) can be solved exactly, with an assumption 
that the major physics would remain unchanged as the 
models are adiabatically tuned to the homogeneous limit. 
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FIG. 1: The sketch of the checkerboard models. The hopping 
signs between dxy and da;z(j,2) orbitals are indicated. 



In this Letter, we present results for three different 
tight-binding models for iron-SCs which are constructed 
using two [28], three [i^, and four orbitals [s^l, respec- 
tively, on a checkerboard lattice (see Fig. [T]). Based on 
the knowledge of exactly solvable four-site problems, we 
show the phase diagrams of the models in the limit of 
weakly coupled plaquettes as a function of intra-orbital 
onsite interaction U , inter-orbital onsite interaction f/'. 
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Hund's coupling and onsite inter-orbital pairing hopping 
Jh: as well as the NNN-AFM coupling J2 which can be 
generated by superexchange mechanism mediated by As 
atoms. In all of the models, the leading SC phase is al- 
ways from the Aig s-wave pairing channel in reasonable 
parameter regions. The superconductivity is intimately 
correlated with the development of the NNN-AFM and 
becomes stronger as J2 increases. Most remarkably, this 
study also shows that the superconductivity and mag- 
netism are orbital-selective in the three- and four-orbital 
models: First, in the three-orbital model, the AFM is 
more pronounced in the d^y orbital while the supercon- 
ductivity is more pronounced in the d^z and dyz orbitals. 
Second, in both three- and four-orbital models, the signs 
of the intra-orbital pairing order parameters of the dxy 
and the d^z (or dyz) orbitals are opposite. The sign dif- 
ference stems from the intrinsic symmetry of the inter- 
orbital hopping amplitudes between dxy and d£cz(i/z) 
This feature could result in the sign-change of the SC or- 
der parameters between the hole Fermi pockets at the F 
point and produce anisotropic or gapless SC gaps in the 
electron pockets around M point in reciprocal space in 
the homogeneous limit. 

Model and method - A generic Hamiltonian of iron-SCs 
can be written as H = Hq + Hj, where Hq is the kinetic 
energy of the d-electrons of irons and Hj describes the 
interactions between them. Explicitly, Hq is given by 

^0 = E E(^^/3 + ^c.S^M.,o.A,P..<r + H-C. (1) 

where a, j3 are orbital indices, i, j label sites and a is spin 
index. Hj can be written as Hi = J2i^io{i) + Hj^, 
where Hio{i) is the onsite interaction, given by. 
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with U, U', Jh, and J{— Jh) denoting intra-orbital 
repulsion, inter-orbital repulsion, ferromagnetic Hund's 
coupling, and inter-orbital pair hopping, respectively. 
H/e describes interactions between different sites. In this 
paper, we specifically consider the NNN (denoted by dou- 
ble angle brackets) magnetic exchange coupling J2 , which 
can be naturally generated by the superexchange mech- 
anism through As atoms. 



Hie — ^ ^ J2^i,a ' Sj_^. 
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Different tight-binding models have been proposed to 
describe the band structures of iron-SCs. In this study, 
we take three different models: a two-orbital model given 
in [2^, a three-orbital model given in [2^, and a four- 
orbital model in [30| . The tight-binding parameters of 




FIG. 2: The pair-binding energy [(a) and (b)] and the diag- 
onal spin-spin correlations [(c) and (d)] for 2x2 plaquette 
in the three-orbital model are plotted as functions of U /W 
and Jh/U with [(a) and (c)] and without [(b) and (d)] ex- 
change coupling J2. The dotted lines in (a) and (b) enclose 
the regions with positive pair-binding energy, suggesting a SC 
phase. 



the models can be found in the above references. We have 
also tried different four-orbital models reduced from the 

All of the models 



five-orbital models constructed in 11 



capture the basic band structures of iron-SCs. In fact, 
the major results reported below are consistent in all of 
these models. 

The checkerboard models are defined on the lattice 
shown in Fig. [1] with the inter-plaquette hopping am- 
plitudes r^''^ and exchange couplings J2 less than intra- 

plaquette parameters, and J2. A four-site plaque- 
tte can be diagonalized exactly, though one may nec- 
essarily resort to the numerics due to its multi-orbital 
complexity. Consequently, we can determine the pair- 
binding energy Ep = 2E{1) - E{0) - E{2), where E{Q) 
is the ground-state energy for a given Q, the number of 
doped holes in an "undoped" reference state. In partic- 
ular, Ep > indicates an effective attraction between 
doped holes. At generic doping, we consider the case 
< {Tjf;,J2} <^ Ep <t: {t^li3,J2}, which allows us to 
obtain a controlled perturbation expansion of the full 
Hamiltonian by small inter-plaquette parameters. Sec- 
ond order perturbation gives us an effective theory de- 
scribing interacting bosons {Q = 2 states) with hopping 
integrals and short-range density-density interactions of 
order /Ep on the effective lattice. Furthermore, it can 
be proved that the ground state enters superfluidity at 
T = 0(31}, except at special doping percentages. Thus, 
the positiveness of Ep could be used as an economic way 
to identify SC phase with generic doping. 
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FIG. 3: The orbital-resolved diagonal spin-spin correlations, 
(Sq -S^) (a, (3 — xz, yz, xy), of a 2 x 2 cluster with 16 electrons 
(2/3 filling) (a) and 14 electrons (b). The parameters are 
given by Jh = 0.2(7, Ji = 0.1(7. Insets: The occupation 
number of different orbitals as a function of U /W with 16 
electrons and 14 electrons, respectively. 



FIG. 4: The pairing amplitude of the intra-orbital (a) and 
the inter-orbital (b) electron pairs as function of U at fixed 
parameters of Jh = 0.2(7 and J2 = 0.1(7. The green dotted 
lines divide the regions of positive (right) and negative (left) 
pairing energy. The 'NN' in the legend means sites 1&2 while 
'NNN' means sites 1&3; all other pairings' amplitude can be 
derived from the Aig symmetry. 



Spin and pairing correlations.- The two-orbital model 
consists of the d^z and dyz orbitals and has been studied 
both analytically and numerically due to its simplicity. 
However, there is no SC phase in the physically mean- 
ingful parameter region in our study if J2 = 0. This 
sharp feature is very similar to the results obtained from 
the functional renormalization group (fRG) for such a 
model[l3|, although the fRG is, in principle, a weak- 
coupling approach. When J2 is added, the supercon- 
ductivity and magnetism in these two orbitals behave 
no qualitative difference from those in three- and four- 
orbital models within the same orbitals. Therefore, we 
will focus on reporting the results from the three- and 
four-orbital models hereafter. 

We summarize the main results for the three-orbital 
model in Figs. [2l |3l and ID Note that the parent com- 
pounds here refer to 2/3 filling, instead of 1/2 filling in 
two- and four-orbital models. In Fig. [21 we plot the pair- 
binding energy and the diagonal (NNN) spin-spin corre- 
lation function as a function of U/W and Jh /U with U' 
satisfying the SU{2) symmetry condition U' = U — 2.5Jh 
and W being the bandwidth. As shown in the figure, 
without J2, there is a small region (circled by the dash 
line) where the SC state is favored and this region is 
located within the region the NNN-AFM is developed. 
With a finite J2 = O.IU, the SC region is significantly 
enlarged as well as the NNN-AFM correlation. It is also 
important to note that a finite Jh is needed in order to 
achieve positive pair-binding energy. 

The NNN spin-spin correlation within a plaquette, 
(g;GS|Sj, • S^|g;GS), as a function of U/W is shown 
in Fig. [3] with fixed Jh = 0-2U and J2 = O.IU. The 
NNN-AFM correlation develops when U reaches close to 
0.4W^ and the intra-orbital correlations are much stronger 
than the inter-orbital correlations. Moreover, the d^y 
orbital has the strongest NNN-AFM correlation as ex- 
pected since the dxy orbital is half filled [see Inset of 
Fig. [nja)]. After doping two holes per plaquette, the 



correlations become even among all the orbitals , as sug- 
gested by equally distributed occupation number among 
all orbitals [see Inset of Fig. [3l^b)]. 

The SC pairing amplitudes, defined as {Q — 
2;GS|A„^(ij)|Q = 0;GS) with A„^(ij) = d,,o.,idj,0,t - 
di,a.fdj,i3^l, are shown in Fig. 131 The positive pair-binding 
energy is achieved when U > Uc ^ 0.4117, close to where 
the NNN-AFM has well developed. As E'p > 0, it is im- 
portant to notice that due to C^v symmetry of a single 
plaquette and the fact |Q = 0; GS), \Q — 2; GS) belong- 
ing to Aig representation, the pair-field operator must 
have precisely Aig s-wave symmetry. The intra-orbital 
NNN pairings in the d^^^^y^'^ orbitals dominate the pair- 
ing strength. The pairing in the dxy orbital is small since 
the dxy orbital remains half-filled and its AFM correla- 
tion is too strong. However, this is an important feature 
that even as the pairing strength is small in the dxy or- 
bital, the sign of the pairing is opposite to the ones in 
the dxz{yz) orbitals. We find that the sign-change is as 
universal as the Aig s-wave pairing symmetry in all of 
the parameter regions we calculate. 

The above major results still hold in the four-orbital 
model as summarized in Fig.jSja). The pairing symmetry 
is still Aig s-wave and with a sign change between the dxy 
and dxz{yz) orbitals. As shown in Fig. |5ja), the positive 
pairing energy is achieved when U reaches about 0.414^ 
as well. The NN-AFM is relatively stronger in the four 
orbital model than the one in the three-orbital model 
since the NNN hopping is weaker. The SC pairing in the 
fourth orbital dx2^y2 is negligible, while the SC pairing 
strength in the dxy orbital becomes stronger, compared 
to the three-orbital model. 

To understand the signs of SC order parameters of dif- 
ferent orbitals, we can consider the effective Josephson 
coupling between them. The Josephson couplings can 
be derived from inter-orbital hoppings. Let's consider 
the inter-orbital hopping between arbitrary two sites (de- 
noted by A and B) described by Hh = ^a/f (^a ^ adB,a,<7 + 
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H.c), where a, /3 are two different orbitals. Assuming it is 
half-filled for both orbitals at the two sites and the spin 
singlet is formed within each orbital due to the intra- 
orbital AFM coupling between two sites, we can easily 
show the inter-orbital hopping can generate the effective 
Josephson coupling between two orbitals given by 



(a) 



H-ff 



j-ABj-AB 



(4) 



where Aq, are the SC order parameters. In the case 
of iron-based superconductors, the general symmetry 
of the lattice requires tt^^ „^ = t!/, ^, and , . = 

—t^,Li„,\ as sketched in Fig. [T] Thus the Josephson 

xz\yz) ,xy 

coupling between the xz and yz is negative while the one 
between the xz{yz) and xy is positive, which favors the 
sign change between A^^^y^^^^^^y^){ij) and A^y^xy{ij). 
This explains the numerical results. 

Discussion.- It is very interesting to compare our 
study to other previous theoretical studies. First, our 
study confirms that the pairing symmetry is dominated 
by the Aig s-wave induced by the NNN-AFM cou- 
pling, which will lead to the sign change between elec- 
tron and hole pockets as many previous studies have 
concluded [3, 

Hill El- Second," our study suggests that 
the Hund's coupling is also important in inducing the SC 
pairing and the NNN-AFM correlation 0. Without the 
Hund's coupling, the model does not lead to the SC in- 
stability. Third, our study indicates the importance of 
the dxy orbital in the strong coupling limit. The pair- 
ing in the d^y orbital favors an opposite sign to those 
in the dyz and dxz orbitals. This result, if still holds 
in the homogeneous limit, may also cause an anisotropic 
SC gap in the electron pockets at the M point of recip- 
rocal space since the electron pockets include both the 
dxy and dxz(yz) orbitals. The presence of dxy can change 
many important properties of the SC pairing. Indeed, 
the importance of dxy orbital has been emphasized in 
the weak coupling approaches, such as the fRG which 
shows that in the absence of the hole Fermi pocket from 
dxy, an anisotropic gap can be developed, similar to our 
implication. However, the physical mechanism, as dis- 
cussed earlier, is completely different from the one in a 
weak coupling approach. Finally, it is worth mentioning 
that the the signs of the Josephson couplings obtained in 
eq|4] explain why the SC state with time reversal symme- 
try (TRS) breaking as suggested in 33, 13] for iron-SCs 
is not favored. A TRS broken state is favored if and only 
if the couplings between all three orbitals are positive. 

Our study provides a new picture of the SC order pa- 
rameters. In the context of iron-SCs, the Aig s-wave 
pairing is generally agreed to exhibit the sign-change be- 
tween electron and hole Fermi pockets via spin fluctua- 
tions. However, our study, though from inhomogeneous 
limit, suggests that in the presence of the dxy orbital, 
there is an additional sign-change between the dxy and 
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FIG. 5: a: The spin correlation (upper) and pairing amplitude 
(lower) of diagonal bonds as functions of U in the four-orbital 
model. The inset shows the pair binding energy in unit of 
band width, b: The schematic pairing amplitude of a 2D 
five-orbital model in the superconducting state inside the first 
(folded) Brillouin zone. Red and blue means positive and 
negative signs, and the thickness indicates the size of the gap. 



the dxz{yz) orbitals. Considering the recent ARPES ex- 
periments, the large hole pocket at the F point is mainly 
from dxy orbital while the small two hole pockets which 
are almost doubly degenerated is composed of the dxz 
and dyz orbitals [35|. Keeping these orbital characters 
in mind, our results do provide a testable, experimental 
prediction: the signs in SC order parameters are opposite 
for the large and small hole pockets at the F point. We 
present a schematic plot for SC order parameters in re- 
ciprocal space in Fig. Mjo)- More physical consequences 
associated with such a sign change will be explored in the 
future study. 
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